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Abstract. We use a dual mesh numerical method to study a non-local parabolic problem 
Ph , arising from the well-known thermistor problem. 
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1 Introduction 

> 

^) \ In this work we propose a dual mesh numerical scheme for analysis of the following non-local 

■ parabolic problem coming from conservation law of electric charges: 

^ : 

\0- ^ - V • {k{u)Vu) = A ^ in Ox]0; Tf 

O; dt Unf{u)dxf 



n = on 9r2x]0;T[, u/t=a = UQ in 17, (1) 

where V denotes the gradient with respect to the x-variables. The nonlinear problem is 
obtained, under some simplificative conditions, by reducing the well-known thermistor problem 
(cf., e.g., |13l 1141 ITC] ^. which consists of the heat equation, with joule heating as a source, and 
subject to current conservation: 

ut = V ■ {k{u)Vu)+a{u)\Vip\^ , V • (cr(n)V(^) = 0, (2) 

where the domain C occupied by the thermistor is a bounded convex polygonal; ip = 
ip{x, t) and u = u{x, t) are, respectively, the distributions of the electric potential and the 
temperature in O; cr(u) and k[u) are, respectively, the temperature-dependant electrical and 
thermal conductivities; (j{u) |V(/3p is the joule heating. The literature on problem ((2j is vast 
(see e.g. 1^1 d El El El El El)- With respect to numerical approximation results to 
problem ((JJ we are aware of ^ IIH I12| I18j : in JHl a numerical analysis of the non-steady 
thermistor problem by a finite element method is discussed; in ^2 the authors study a spatially 
and completely discrete finite element model; in a semi-discretization by the backward Euler 
scheme is given for the special case k = Id] in ^ a box approximation scheme is presented and 
analyzed. A completely discrete scheme based on the backward Euler method with semi-implicit 
linearization to ^ is presented in ^21 for the special case k{u) = 1. Existence and uniqueness 
of solutions to the problem Q were proved in ^U] . 
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Finite volume methods emerged recently and seem to have a significant role on concrete 
applications, because they have very interesting properties in view of the subjacent physical 
problems: in particular in conservation of flows. An equation coming from a conservation law 
has a good chance to be correctly discretized by the finite volume method. We also recall that 
these schemes have been widely used to approximate solutions of the heat linear equation, semi- 
linear or parabolic equations. Since we consider data / with lack of regularity when compared to 
previous work, we need a new way to discretize We present a dual mesh method capable of 
handling the non local term - — '^■^^"^ which is a noticeable feature of d), by generalizing the 

results of A box approximation scheme for discretizing with the case k being different 
from the identity is obtained. Speed of convergence is directly related with regularity of the 
continuous problem. When one increases regularity of the second term and data, the solution 
see its regularity increasing in parallel, and precise speed of convergence can be established. 
In the existing literature (see e.g. El) the error estimates for both the finite element or 
volume element method are usually derived for solutions that are sufficiently smooth. Because 
the domain is polygonal, special attention has to be paid to regularity of the exact solution. 
We give sufficient conditions in terms of data and the solution u that yield error estimates (see 
hypothesis (HI) below). 

The text is organized as follows. In Section El we set up the notation and the functional 
spaces used throughout the paper. Section |21 introduces a box scheme model for problem 
and existence and uniqueness of the solution of the approximating problem (|12j) is obtained from 
the fixed point theorem and equivalence of norms in the finite dimensional space Sf^. Finally, 
in Section m under some regularity assumptions, we prove error estimates. 



2 Notation and functional spaces 

Let (•, •) and ||-|| denote the inner product and norm in L^(J7); Hq{^) = |ii € H^{^), u/dQ. = O}; 
II • lis, II • \\s^p denote the H^{Q) and the W'^'^{Vl) norm respectively; Th denote a triangulation 
of Q.; be the set of vertices of a quasi-uniform triangulation T^; and {'S'{]}^^q be the family 
of approximating subspaces of Hq{^) defined by 

Sj] = € Hq{Q.) : u/e is a linear function for all e G T/i} . 

In the remainder of this paper we denote by c various constants that may depend on the data 
of the problem, and that are not necessarily the same at each occurrence. We assume that the 
family of triangulations is such that the following estimates Hj hold for all f G SjJ: 

lbll/3,g < c/i""'^~2'^^"^°'^/P~^/^^||7;||,,p, 0<r </?<!, l<p,g<cx), 

lb ||o,oo < c| ln/l|5||i;||i. (3) 

Let P/j : L^(r2) — > Sf^ be the standard L^-projection. One has j^: 

lit; - Phv\\ + h\\v - Phv\\i < ch^\\v\\2, 

\\v - Phv\\o,oo < ch\\v\\2, ll-Ph-vlli.oo < c||i;||i,oo- (4) 

We construct the box scheme Bf^ (dual mesh) employed in the discretization as follows. From 
a given triangle e G T/^, we choose a point g' G e as the intersection of the perpendicular bisectors 
of the three edges of e. Then, we connect q by straight-line segments to the edge midpoints of e. 
To each vertex p G T^, we associate the box bp G Bh, consisting of the union of subregions which 
have p as a corner (see Fig. For the piecewise constant interpolation operator I^, defined by 

Ih:C{n)^L\n), IhV = v{p), on bp G B^, ^peT^, 
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Figure 1. Construction of the dual mesh. 



we have the following standard error estimates ^ : 



\v\\ < \\Ihv\\ < c\\v\ 



\v — Ihv\\ < c/i||f 111, V -y G 5| 



• 



(5) 



We denote by Nh{p) the set of the neighboring vertices of p G T^, db = UpGT'' ^^P' ~ 
\-\p*£Nh,{p)i^pp*}' "^here Tpp* = dbpf]dbp^, (see Fig. [T]). Let Igi, : 56 ^ ]R+ be defined as follows: 
for p & and bp € Bh, 

Idb/^pp* = \P-P*\ for p* G Nh{p). 

For b G Bh, we denote the jump in w across db at x by [w\Qh{x) = w{x + 0) — w{x — 0), where 
w{x lb 0) are the outside and inside limit values of w{x) along the normal directions for db. 

We now collect from the literature some important lemmas and trace results, that are 

needed in the sequel. 

Lemma 1. Assume that Bh is a dual mesh. If v is a piecewise linear function, and x is not 
a vertex, then 



[hv]/dbj,{x) 



dv 



Idb/^f 



X G 



ybeBh, 



where n is the unit outward normal vector on db. 
The /i-dependent norms are defined as follows: 

and \\v\\o h = WhvW ■ 



V hh 



\ledb 



Lemma 2. There exists a constant c > such that 

c^^ \\vv\\ < \\v\\i^h < c\\vv\\, y V e Sh, 

c~"^||f|| < ||^^||o,/i < cllwll, V w G S^. 
Lemma 3. For any a G C{^}) there exists a positive constant c such that 

< c||ti||i||f||i, Vn,t>G5/]. 



E 

beBh 



du 
db dn 



Moreover, if there exists a constant oq > such that a > oq inVl, then 

dv 



<;-'llt.|l?< 



E 



a^hv, y v £ Sh. 
db dn 



(6) 



(7) 
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Let Qh '■ H'^{^) — > S'^ be defined by QhU — ihU € S^, and 



^ / d{u- Qhu) . n w ^ cO /e^ 



where i/i : C{VL) is the Lagrangian interpolation operator and u G H'^{Q). 

Lemma 4. Assume that a S wii/i a > oq /or some constant oq > 0. T/ien, i/iere exisis 

c > stic/i i/iai /or u e H'^{^}) 

\\u - Qhu\\i < ch\\u\\2. (9) 
Moreover, i/ u € p| p^/i.oo^j^)^ ^/^g^ 

IIQ/i-wlli.oo < c(||n||i,oo + \\u\\2) ■ (10) 
Lemma 5. For each b G -B/j one has 

Throughout this work, we assume that the following hypotheses on the solution and data of 
problem are satisfied: 

(HI) u € L^{Hl{n){^H^{^)), ut e L\H\n))- 
(H2) < k{s) < c; 

(H3) there exist positive constants ci, C2 and z^, such that u < /(^) < ci|^| + C2 for all ^ G M; 
(H4) \f(0-f(C)\ + \m-kia<c\C-a 

3 Existence and uniqueness result for the box scheme method 

Let u be the solution of Integrating over an element b in S/^ we obtain: 

[u,- [ k{u)^= \ J fiu), ybeB,. (11) 

Jb Jdb on f{u)dx) Jb 

We consider a box scheme defined as follows: find G Sf^ such that 

(4»f.4.)-E//MlS,^--(j;7^(/M.M. (12, 

where m'*(0) = PhUo and Ih is the interpolation operator. 

Theorem 1. Let (H1)-(H4) be satisfied. Then, for each h > 0, there exists tQ{h) such that 
possesses a unique solution for < t < to^h). 

Proof. We begin by proving existence of solution. We define a nonlinear operator G from S'jl 
to as follows. For each G S^, = G{u^) is obtained as the unique solution of the 
following problem: 



beBh 
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We remark that G is well defined. Using v = as a test function in H13|) . hypotheses (H2) 
and (H3), and Holder's inequality, we can write: 



l^Wl^yj^^ + cWvjX < C{f{u^),hw'^) <cj {\U'^\ + 1)14 

< c\\u''\\L2\\Ihw''\\L2 + cWhw'^'W < c||n^||i||4«;^||i + c\\hw'' 



<^\\I^ro^\\l + c\\u^\\l + c. 
Thus, we have 

+ c||«;'^||?<c||n'^ll? + c. (14) 

Integrating ((TUl with respect to t and using the equivalency of \\Ih • \\ and || • || in S'^ (see ((SJ) 
yields 



+ c [ \\w''\\l < c\\IhPhUof + c [ \\u''\\ldx + ct 
Jo Jo 

< c||tio||^ + c / + ct. 







Define now the following set 

D=L''e Si Wu^'f + c I \\u^\\j< c{\\uor + 1] 



We can easily see that D is closed subset of L°°(0, t, L^(r2)) with its natural norm. We conclude 
that there exists t > such that G{D) C D. To obtain that G has a fixed point = G{w^), 
we prove that G is a contraction. Conclusion follows from Banach's fixed point theorem. For 
this purpose, let and € 5*° x 5° such that Gu^ = and Gu2 = W2- We have, from the 
equation (fT3|) verified by w\ and that 

fetX Jdb Jdb dn 

^ (/ (4) , 4^) - {f{u'i)jHv). 



On the other hand, one has 



Schwartz inequality implies that 



Wl - U2\\o,db I lit"! 



baBu-""" \beBh 

By Lemma El we have 

beBh b£Bh 
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Thus, from the inverse estimate (jSJ 
On the basis of hypotheses (H1)-(H4), we have: 



j^^ [KO - H^2))^hv > -cmw%\\ui - u%\\v\\,. (15) 



^ {f{u^,),hv)- „ (/ (4) , hv] 



^ 2ifi^l)-fi^2),hv) 



{In f i^i) dx) 



(L / (^2) d^) Un f iO dx) 

< c\\u'l - 4|| \\v\\ + c||4w||i2(n)||'Ui - u^Wmn) 

< c\\u'l - u^W \\v\\ < c||uj - u^lli (16) 
It follows from and that 

(4«-u;,^),/..)- / 

<c(/i)||<-4||i||Hli- E / {H^i)-Hu!^))^hv<c{h)\\u'l-u%\\v\\i. (17) 

Now, using V = w'l — W2 as a test function in 1)17^ . we obtain from ((T)): 

- ^2)\\' + c\H - w^Wl < cih)\\u1 - u%\\w'l - w%. (18) 

With use of the Holder's inequality and equivalency of • || and || • ||, integration of ((TH|) with 
respect to time gives: 



\{w1-w!^)f <c\\h{w1-w!^)f <c{h) [ \\{4-4)\\lds 



Thus G is a contraction. We prove now uniqueness. Following the same arguments as before, 
we have 

{lh{4t-ni),Ihv)-Yl I fc(n?)^^=^4i;<c(/i)||4-4llilbll- (19) 
Choosing v = Ui — U2 as test function in H19|) . using again Q and integrating, we obtain 

||K-4)||'<c(/i) r||(n?-4)||'d5, 
Jo 

which gives, by Gronwall's Lemma, uniqueness of solution. ■ 
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4 Error analysis 



In this section we prove error estimates under certain assumptions on regularity of the exact 
solution u. 

Theorem 2. Under assumptions (H1)-(H4), if (n, n'^) are solutions of ()11() - ()12() for <t < 
to{h), then 

Proof. From and we obtain 

fees. 



dn 



X 



X 



-Y. f {Ky)-k{u''))^hv+{{I-Ph)uuhv) + {{I-Ih)PhUuhv)- (20) 

We now estimate, separately, the terms on the right-hand side of (|2fl|) . We have from © and Q) 
that 



b&Bh 



db 



dn 



< c\\PhU — < c/i||ii||2 ||'y||i < c/i||f 111, 



u — u 



du 



dn 



) I <c/ii||^||i ( 



\u-u^\\ldb 



By Lemma 13 inverse inequality ^ and @, we have: 

11^'^ - <,db < 2 ||n^ - P^n||2 + 2 ^ W^hU - u\\Iq, 

beBh beBh b£Bh 

<c{h^ + \\u'' -Phuf). 



(21) 



(22) 



(23) 



Consequently, we obtain from (|23() that 



< c(/i+ ||n^ -P/,n||)||7;||i. 



Based on our earlier development in ()16p . we also know: 



\f{n''),hv) 



{f{u),Ihv) < c\\u - n||||f 111. 



(24) 



{f^f{u^)dxy' {f^f{u)dx) 

Let V = — PhU be a test function in 1)201) . Using Lemma 01 it follows from (|21() -()25 p that 
1 d 



(25) 



2dt 



Whiv!" - Phu)\\^ + cWu'' - PhuWl 
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< c\\u'' - u\\\\u^ - Phu\\i + ch\\u^ - Phu\\i + c{h + \\u^ - Phu\\) \\u^ - Phu\\i 
+ c(||(/ - Ph)ut\\ + 11(1 - h)Phnt\\)\\u'' - Phu\\ 

< dWil - Ph)Ut\\ + 11(1 - h)PhUt\\)\\u'' - PhU\\ 

+ c{h + ||n^ - Phu\\) \\u^ - PhuWi + \\PhU - u\\\\u^ - Phu\\i. 

By properties ^ and Cauchy's inequality, it follows: 

- Phu) f + c\\u^ - PhuWl < c\\PhU - ^x|li||n^ - P;,n||i 
+ c{h + 11(7 - Ph)ut\\ + 11(7 - Ih)PhUt\\ + - Phu\\}\\u'' - PhuWi 

< c{h^ + 11(7 - Ph)ut\\l + 11(7 - h)PhUt\\^ + Ik'^ - Phuf} + - PhuWl 

< c{h^ + h^WutWl + ch^WPhUtWl} + c||n^' - P,nf + ^Wu"^ - P,n||2. 

Hence, 

j^\\h{n'' - Phu) f + c\\u'' - Phu\\\ < ch^ + c\\u'' - Phuf. (26) 

Integrating (|26|) and applying Gronwall Lemma and using again the equivalency of || • || and 
\\Ih • II, we get that 



- Phu\\^ + c [ ||n^ - PhuWl < ch\ 
Jo 



Then, by the triangular inequality, we conclude with the intended result. ■ 

Under more restrictive hypotheses on the data, it is possible to derive the following error 
estimate. 

Theorem 3. Assume (H1)-(H4). If k{s) = 1 and uq € H^{n)f]H'^{n), then 

~ ^IIl°°(_H'1) ^ c/l. 
Proof. From equations 1)11(1 and ()12() . we have: 

{hu1-ut,Ihv)-Yl [ ^hv+Y. I ^^^"^ 

[kfi'^ndx) [\^f{u)dx) 
Using the definition (jH)) of Qh, we get 

d {u^ - Qhu) 



{hn, -nuhv) - Y ^ hv 



^ {f{u''),hv) - -—^--^{f{u),hv), 



and it follows that 

bdBh''^^ 
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+ ((/ - Qh)ut,Ihv) + ((/ - Ih)QhUt, hv). 



(27) 



In order to estimate the right hand side of the last inequality, we treat both terms separately. 
By similar arguments to those used in (|T6]) . 



A 



A 



{f{u)Jhv] 



< cWv!^ — u\\ \\v\ 



Taking a function test v = {u^ — Qhu)t in by Q and (fTn|) we have 



\\Hu^-Qhu)tf- I 



d{u^ - Qhu) ^ f h n \ 
h{u -Qhu)^ 



(28) 



< c{h lh)QhUt\\ + Qh)ut\\ + WQhu - u%\\h{u^ - Qhu)t\\ 

< c{h + chWuth + chWQhUtWi + \\QhU - u''\\}\\Ih{u'' - Qhu)t\\ 

< c{h + WQhU - u''\\}\\Ih{u^ - Qhu)t\\ 

< ch^ + c\\QhU - u'^f + ^||//^(u'^ - Qhn)tf. 
Integrating ()28p . we arrive to 

<c{h^ + \\u^ - Qhuwi^ = ch^ + c\\u" - QhuWli^H^Q))^ 
and Theorem [21 gives 

\\u^-Qhu\\l < ch^- 

On the other hand, by triangular inequality, Q and the regularity of the exact solution u, we 
have 

\\u^ - u\\l < 2\\u'' - Qhu\\l + 2\\QhU - u\\l < ch^\\u\\l + ch^ < ch^. 
We conclude then with the desired error estimate. ■ 



5 Conclusion 

In this paper a dual mesh numerical scheme was proposed for a nonlocal thermistor problem. 
We have showed the existence and uniqueness of the approximate solution via Banach's fixed 
point theorem. We have also proved iJ^-error bounds under minimal regularity assumptions. 
We only obtain first-order estimates: higher order estimates are difficult to obtain due to the 
nonstandard nonlocal term. Optimal error analysis to the present context, under appropriate 
smoothness assumptions on data, can be derived by application of the techniques of [H] , but this 
needs further developments. 
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